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Contributions of three-nucleon forces (3NF) to proton-deuteron scattering observables at energies 
below the deuteron breakup threshold are studied by solving the Faddeev equation that includes the 
Coulomb interaction. At E p = 3.0 MeV, we find that the central part of a two-pion exchange 3NF 
removes the discrepancy between measured cross sections and the calculated ones by two-nucleon 
forces, and improves the agreement with T22 experimental data. However, the tensor part of the 
3NF fails in reproducing data of the analyzing power T21 by giving worse agreement between the 
measured and the calculated. Detailed examinations of scattering amplitudes suggest that a P-wave 
contribution in spin quartet tensor amplitudes has unsuitable sign for reproducing the T21 data. 

PACS numbers: 25.10.+S, 21.30.-x, 24.70.+S 



Contributions of three-nucleon forces (3NF) have been 
studied extensively for proton-deuteron (pd) scattering 
observables, since the interactions are successful |l| in 
solving the problem of three-nucleon (3N) underbinding 
for realistic two-nucleon forces (2NF). However, the Fad- 
deev calculation, which is one of the practical methods 
to treat 3N systems, has conventionally neglected the 
Coulomb interaction because of mathematical difficulties, 
while the interaction is essentially important in the low- 
energy pd scattering. Recently the problem has been 
solved by an Faddeev integral equation approach at en- 
ergies below the deuteron breakup threshold @], where 
the phase-shift parameters by the calculation agree to 
those by an Faddeev differential equation approach Q 
as well as the Kohn variational method 0, 0| with very 
good accuracy. This allows us to investigate the con- 
tribution of 3NF in the low-energy pd scattering by the 
Faddeev calculation including the Coulomb interaction. 
In this Rapid Communication, we will report the main 
result of such calculations, where one will find remarkable 
effects of the 3NF with promising success of the central 
part but pessimistic failure of the tensor one. 

In the calculation, we adopt the Argonne Vis model 
(AV18) @ for the 2NF and the Brazil model (BR) for 
the two-pion exchange 3NF. Further we introduce two 
kinds of 3NF: a spin-independent Gaussian (GS) 3NF Q 
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where p 2 = §(rf 2 + r| 3 + r%), a = 1.5 fm -1 , W = -20 
MeV, and Pn is the projection operator to the spin and 
isospin triplet state of the pair. Calculated 3 He 

binding energies are 7.79 MeV for the BR-3NF in addi- 
tion to the AV18 (AV18+BR), 7.74 MeV for the GS-3NF 
(AV18+GS), and 7.74 MeV for the BR- and SO-3NFs 
(AV18+BR+SO), which are compared with the empiri- 
cal value of 7.72 MeV. The GS-3NF simulates the central 
part of the BR-3NF @, and then the difference between 
both calculations describes the contribution of the ten- 
sor part of the BR-3NF. The SO-3NF is adopted as a 
simulation of the spin-vector type 3NF that reproduces 
empirical vector analyzing powers, but the origin is yet 
unknown at present. 

The comparison between the calculated quantities and 
the measured ones [RJ, [ll|, [l2| is shown in Figs. ^ 121 an d 
Elfor E p = 0.65, 2.5, and 3.0 MeV, respectively, where 
the differential cross section da/dQ, the vector analyz- 
ing power of the proton A y , that of the deuteron iTn, 
and the tensor analyzing powers T20, T21, and T22 are 
displayed. Although overall agreements are obtained be- 
tween the calculations and the experimental data at the 
three incident energies, the BR-3NF seems to deterio- 
rate the agreement to the T%\ data of middle angles at 
E p = 2.5 and 3.0 MeV. This effect can be confirmed by 
X 2 /datum in Table[IJ where we include similar analysis of 
the data [13j at E p — 1.0 MeV to see energy dependence 
of x 2 in detail. 

To demonstrate characteristic features of such BR-3NF 
contributions, we will show, in Fig. ^ da/dfl, T20, T21, 
and T 2 2 typically at E p = 3.0 MeV, divided by the theo- 
retical values obtained by the AV18 calculations. In the 
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FIG. 1: Differential cross sections da/dQ., vector analyzing 
powers A y and iTu, and tensor analyzing powers T20, T21, 
and T22 of the pd scattering at E p = 0.65 MeV. The data are 
taken from Ref. [Toll . The solid, dashed, dotted, and dash- 
dotted lines are the calculations by the AV18, the AV18+BR, 
the AV18+GS, and the AV18+BR+SO, respectively. 



figure, the solid horizontal lines (the 2NF-lines) describe 
the AV18 calculations, and the deviations of the theoret- 
ical curves from the 2NF-lines describe the 3NF contri- 
butions. The experimental points of da/dfl deviate from 
the 2NF-line with characteristic angular distribution, and 
the deviation is well reproduced by the calculations in- 
cluding the 3NF effects. This 3NF effect is attributed to 
the central part of the BR-3NF, because the calculations 
by the AV18+BR and those by the AV18+GS give sim- 
ilar results, and the additional SO-3NF produces very 
small contributions. As was discussed in Ref. [2|, this 
3NF contribution is produced through the spin doublet 
scalar amplitude and is related to the 3NF contribution 
in the 3N binding energy [l4j ]. 

Contrary to such success, the BR-3NF fails in repro- 
ducing the experimental data of T 21 . Fig. 21 shows that 
T 21 calculated by the AV18+BR at middle angles is lo- 
cated in the opposite side of the measured values with 



FIG. 2: Same as Fig.[T]but at E v = 2.5 MeV. The data are 
taken from Refs. [TllllSj. 

respect to the 2NF-line, indicating the opposite sign of 
the 3NF contribution to be desirable for reproducing the 
data. This failure is dominantly due to the tensor part of 
the BR-3NF, since the contributions of the GS-3NF and 
the SO-3NF are small. Similar discrepancy is also ob- 
served in T20, but the role of the tensor part is not clear. 
The agreement with the T22 data is improved, though 
is not sufficient, by taking account of the 3NF, where 
the dominant contribution of the 3NF comes from the 
central part of the BR-3NF since three calculations by 
the AV18+BR, the AV18+GS, and the AV18+BR+SO 
calculations give similar results. 

Next we will examine the 3NF contributions in more 
detail for T21 and T22. To understand roles of the central 
interaction and the tensor one of the 3NF individually, 
we will decompose the scattering amplitude according to 
the tensor property in the spin-space by expanding the 
T-matrix M into the spin-space tensors Si 

M = j2(-rs {K K ) Ri K \ (3) 

Kk 

where Rk is the coordinate-space tensor and K (k) is 
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FIG. 3: Same as Fig. [Tj but at E p = 3.0 MeV. The data are 
taken from Refs. [nfll^. 



the rank {z component) of the tensor. Then the scatter- 
ing amplitude is given by Ref. |l5l | as 

(v' p v' d ;k f \M\u p v d ;ki) 
= ^{spSdV p Vd\s l v l ){s p s d v' p v' d \s } v } ){-) s f- v f 

s iSf 

xY,{^ fVl - Vf \KK)M^\ Sl s f ) (4) 
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where K = K (K + 1) for K = even (odd), and Q m (fc) 
is related to Yg m (k) as usual 0- The quantity Sj (s/) de- 
notes the channel spin in the initial (final) state, which 
is \ (the spin doublet states) or | (the spin quartet 
states). The function F^ K \siSf£i), the invariant am- 
plitude, is a function of the scattering angle 9 and is 



TABLE I: X 2 /datum of the AV18, AV18+GS, AV18+BR, 
and AV18+BR+SO pd observables compared with the exper- 
imental data from Refs. [TollTl ITllllt . 

a A y iTn T20 T21 T22 



£ p =0.65 MeV 
AV18 21 16 

AV18+GS 1.7 12 
AV18+BR 2.0 12 
AV18+BR+SO 1.7 2.4 



13 2.9 3.4 2.7 

7.8 2.3 2.7 2.0 

9.5 2.6 3.1 2.1 

1.5 2.6 3.2 2.1 



E p =1.0 MeV 

AV18 51 164 60 

AV18+GS 2.8 115 28 

AV18+BR 2.8 116 37 

AV18+BR+SO 2.5 20 2.1 



3.7 3.8 4.2 
1.4 2.8 1.5 
1.6 3.4 

1.8 3.5 



1.7 
1.8 



AV18 23 189 90 1.5 5.2 10 

AV18+GS 1.9 111 33 3.7 5.0 2.6 

AV18+BR 1.8 111 53 6.3 10 2.8 

AV18+BR+SQ 2.2 1.9 10 5.5 8.9 3.1 
E p =3.0 MeV 

AV18 25 160 114 3.2 7.1 18.0 

AV18+GS 2.3 95 47 5.1 8.3 3.6 

AV18+BR 2.2 94 74 10 20 4.3 

AV18+BR+SO 2.6 1.9 6.0 9.0 18 5.3 



designated by the tensor rank K. Thus the amplitude 
describes the scattering by interactions classified by K, 
i.e., F^°\siSf£i) {F^ 2 \siSfli)) describes the scattering 
by the central (tensor) interactions. Because of the time 
reversal theorem 

F {K) {s fSi l f ) = (-r- s fF^(s iSf k), (6) 

only five amplitudes of nine tensor ones are independent 
0. 

In the low-energy scattering, the scalar amplitudes 
dominate other amplitudes as shown in Fig. where 
the magnitudes of the two scalar and five independent 
tensor amplitudes at 8 = 90° are displayed as functions 
of E p . Thus we will apply to the tensor analyzing powers 
T2k{k — 0, 1, 2) an approximation in which terms not in- 
cluding the scalar amplitudes are neglected, and get the 
following expressions: 
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FIG. 4: Differential cross sections da/d£l and tensor analyzing powers T20, T21, and T22 at E p = 3.0 MeV divided by the 
theoretical values with the AVf8. See the caption of Fig.Qfor the definitions of the theoretical curves. 



with 

N R = Tr(MMf) - 6^. (11) 

Now our attention will be focused on the special scat- 
tering angle 9 = 90° for simple considerations. Each 
component and the sum in Eq. J7J) at 9 = 90° are dis- 
played for T21 (T 22 ) in Fig. E^a) (Fig. HJb)) as func- 
tions of E p . Fig. EJa) shows that the unfavorable ef- 
fect of the BR-3NF on T21 seen in Fig. 0] appears for 
E p > 2 MeV, and the largest contribution arises from 

[31 

T^ J of the spin quartet scattering. Since the quartet 

scalar amplitude Mj°^(||) is hardly affected by the 3NFs 
as shown in Ref. the spin quartet tensor amplitude 

m| 2 -*(||) included in should be responsible for the 
effect. Adopting the Madison convention for the coordi- 
nate system, z \\ ki and y || fej x kf, one can write down 

A ( 2) 

[C^(fci) (g> Ci f (kf))\ ' as functions of cos0 and sin#. Re- 
ferring to Eq. JSJ), we can restrict the invariant amplitude 
effective for M| (sjS/) to ii = if = 1 (P-wave) as long 
as we are concerned with 9 = 90°, since the amplitudes 
of (£i,£f) = (0,2) vanish because of the factor cosflsin^, 
and those of (£i,£f) = (2,0) do not appear. In Eq. JSJl, 
ii {If) gives the orbital angular momentum in the initial 
(final) channel, when the ^-dependence of F^ K \siSf£i) 
is neglected because of low-energy scattering. Then we 
conclude that the P-wave scattering in the spin quartet 
state has a key of the unfavorable 3NF contribution to 



T 2 i. 

The 3NF contribution to T22 will be analyzed in a sim- 
ilar way. As seen in Fig. EJb) , gives small contribu- 

[21 [31 

tions to T22. Other two terms T 22 and T 22 are accompa- 
nied by the scalar amplitude Mq°'(||) as Eqs. © and 
HI U| l. which receives effectively no 3NF contribution as 
discussed above. The related tensor amplitude M^\~) 
and M 2 2 ^(||) are also scarcely affected by the 3NF. In 

fact, M 2 2 ^(||) has been shown to be unaffected by the 
3NF in neutron-deuteron scattering at E n — 3 MeV 
However, Nr in the denominator is affected by the 3NF, 
and most of 3NF contributions to T22 are produced by 
this effect. Since Nr is proportional to the cross sec- 
tion, the successful improvements by the 3NF in da/dfl 
and T22 are achieved due to the same origin, the central 
component of the 3NF. 

We will conclude that at E p = 3.0 MeV, the central 
part of the BR-3NF produces the successful contribution 
on the differential cross section as well as on the T22 an- 
alyzing power, while the tensor part of the interaction 
gives the undesirable contribution to T21. The 3NF ten- 
sor effect on T21 is not observed below 1 MeV as seen in 
Table U and Fig. El suggests it may become appreciable 
above 2 MeV. Thus precise measurements of the observ- 
ables for the pd scattering are highly desirable in such an 
energy region for the study of the 3NF tensor effect. Also 
it will be interesting to examine the energy dependence 
of the effect up to higher energies, where several prob- 



5 



lems of the tensor analyzing powers have been reported 
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FIG. 6: Contributions of (90°), (90° ) , and T 2 ( ^(90°) 
to T2 K (90°) as functions of E p for k = 1 (a) and k = 2 (b). 
The dashed lines, the solid ones, and the square blocks rep- 
resent the calculations by the AV18, the AV18+BR, and the 
AV18+GS, respectively. 



